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Red-giant stars are emerging as one of the most interesting areas of space asteroseismology. 
Even a relatively basic analysis leads to the determination of the global parameters of the stars, 
such as their mass and radius, and the very extensive space-based data now available for a 
large number of stars allow detailed investigation of the deep interiors of red giants, including 
distinguishing between stars that do and do not have helium fusion in the core, on the basis of 
periods of gravity waves partially trapped in the core. Here I review the theoretical background 
for these new developments and provide a simple explanation for the effect on the period spacing 
of central helium fusion. 



1. Introduction 

Asteroseismology, and hence the study of stellar propertie s, is being revoluti onized by 
the extremely accurate and extensive data from the CoRoT (jBaglin et al. I l2009l) and Ke- 
pler (Borucki ct al. 2009) space missions. Analysis of timeseries of unprecedented extent, 
continuity and sensitivity have allowed the study of a broad ran ge of stellar variabil- 



ity, including oscillati ons of a v ariety of pulsating stars (see, e.g.. iGilliland et al. 2010l : 
Christensen-Da lsgaa rd fc Thom pson 201 ll, for reviews), and leading to comparative as- 
teroseism ology (or synasteros eismoloq'^ for main-sequence stars showing solar-like os- 
cillation ([Chaplin et al.ll2011l ). Also, early analyses of Kepler data have demonstrated 
the power of asteroseismology in characterizing the ce ntral stars in planetary systems 
( Christensen-Dalsgaard et al. 2010l : Batalha et al.ll2011 ) , and such investigations will un- 
doubtedly play a major role in the continuing Kepler exploration of extra-solar planetary 
systems. 

However, perhaps the most striking results of space asteroseismology have come from 
the investigation of red giants. Given the extensiv e outer convection zones of red gi- 
ants, solar-like oscillations were predicted quite early (jChristensen-Dalsgaard fc Frandsen _ 
1983 ) . Ground-based observa tions have been carried out in a few cases (e.g.. lFrandsen et al, 
20021: iDe Ridder et al.|[2006l) but owing to the very long periods of these huge stars such 
observations are extremely demanding in terms of observing and observer's time. Space 
observations, on the other hand, allow nearly continuous observa tions over very exte nded 
periods, as demonstrated b y early observations by the WIRE ( Retter et al. 20031 ) and 
MOST (jBarban et al.ll2007h ' satellites. With CoRoT and Kepler we have obtained exten- 
sive data for literally thousands of red giants, providing precise characterization of their 
overall characteristics, as well as, in some cases, detailed information on the properties 
of their deep interiors. This promises completely to change our understanding of these 
late phases of stellar evolution. 

Here I provide an overview of the oscillation properties of red giants and discuss some 
of the initial results of the analysis of data from CoRoT and Kepler. An introduction 
to the properties of solar-like oscillations, emphasizing the observational aspects, was 
provided by Beddi ng (this volume). Furth er introductions to solar-like stellar oscillations 



were provided by IChristensen-Dalsgaard (|2004) and, in particular, by the monograph 



f a term coined by D. O. Gough 
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by lAerts et al. (2Ql3)- To illustrate the relevant properties, computations of red-giant 
models and oscillations are presented. These probably do not fully cover the complex- 
ities of these late stages of stellar evolution and hence are not intended as models of 
specific observations, but it is hoped that they at least capture the main features of these 
stars. The modelling was carri ed out using the Aarhus STellar Evolution Code AS TEC 
( Christensen-Dalsgaard 2008a), w ith OPAL equat ion of state ( Rogers et al. 19961) and 
opacities (jlglesias fc Rogers! Il996l) . and NACRE jAn gulo et alj [l999|) nuclear reaction 
rates. Convection was treated with the iBohm-Vitense ( 1958 ) mixing- length formalism, 
with a mixing length proportional to the pressure scale height and approximately cal- 
ibrated to a solar model. Diffusion and settling were neglected and no mass loss was 
included. Adia batic oscillations were calcula tion using the Aarhus adiabatic pulsation 
code ADIPLS ( Christensen-Dalseaar dl l2008bh . To resolve the very rapid variations in the 
eigenfunctions in the stellar core, the models were transferred from the original mesh 
used in the model calculation to a much finer mesh, designed to capture adequately the 
properties of the eigenfunctions. However, it is probably fair to say that the mesh distri- 
bution, and more generally the numerical procedures to deal with this type of oscillation, 
requires further optimization and testing. 



2. Simple properties of stellar oscillations 

This section provide a brief background on the theory of stellar oscillations, as far as 
the oscillation frequencies and their relation to stellar properties are concerned. Some 
emphasis is given to the oscillation properties of red giants. 

To understand the properties of stellar oscillation frequencies, it is a good approxima- 
tion to assume that the oscillations are adiabatic, and hence ignore the processes that 
damp or excite the modes. This approximation breaks down in the near-surface layers, 
where there are strong departures from adiabaticity. Also, the modelling of this region 
is highly uncertain, owing to the detailed effects of convection on the mean structure of 
the star and the dynamics of the oscillations. These near-surface processes have a sig- 
nificant effect on the frequencies, as is evident in the analysis of solar observations, and 
they must be tak en into account in de tailed asteroseismic analyses of solar-like oscillation 
frequencies (e.g., Kieldsen et al. 20081) . However, for the present discussion we can ignore 



them. 

With this approximation the problem of computing stellar oscillation frequencies for a 
given stellar model is relatively straightforward Q This feature, together with the ability 
to determine the observed frequencies with very high accuracy, makes stellar oscillations 
such excellent diagnostics of the stellar interior. For the purpose of computing adiabatic 
frequencies, and assuming that the equilibrium model satisfies hydrostatic equilibrium 
(see Kawaler, this volume), the structure of the star is characterized by the density p{r) 
as a function of distance r to the centre, as well as the adiabatic exponent Ti(r), where 

ri-fS^l , (2.1) 



dlnp 



ad 



where p is pressure and the derivative corresponds to an adiabatic change. Given p{r) 
the mass distribution in the star can immediately be determined, and p{r) can then 
be computed from the equation of hydrostatic equilibrium, with a suitable boundary 
condition. This essentially defines all that is needed to complete the equations of adiabatic 
oscillation. For the predominantly acoustic solar-like oscillations the most important 



f Although, as discussed below some care is needed when considering red giants. 
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quantity is the adiabatic sound speed c, given by 

2 Tip 



(2.2) 



In many relevant cases, the equation of state can be approximated by the ideal gas law: 



P 



(2.3) 



where T is temperature, fj. the mean molecular weight, Boltzmann's constant and niu 
the atomic mass unit. We then obtain 



rifceT 



(2.4) 



We consider small-amplitude oscillations in a slowly rotating star. Then the dependence 
of the eigenfunctions on co-latitude 9 and longitude (/) can be separated as spherical 



harmonics Y,^' 



the degree I measures the total number of nodal lines on the stellar 



surface and provides a measure of the local horizontal wavenumber /ch^ 



Kjj — 



1(1 + 1) 



(2.5) 



while the azimuthal order m, with |m| ^ /, measures the number of nodal lines crossing 
the equator. For each (l.,in) the star has a set of eigenfrequencies uJnim labelled by the 
radial order n. For adiabatic oscillations ujnim is real; the oscillation depends on time t and 
(p as cos{m(p — ujt) , i.e., for m 7^ 0, as a running wave in longitude. For a nonrotating star 
the frequencies are independent of m. Slow rotation with angular velocity introduces 
a frequency splitting 

— ^nlO + rnPnlm{^)nlm , (2.6) 

where {Q)nim is an average over the star with a weight determined by the oscillation 
eigenfunction. If w = ^(r) depends on r alone finim = Pni and the average are indepen- 
dent of m and the rotational splitting is linear in m. 

After separation of 1^'" the equations of adiabatic oscillations reduce to a set of four 
linear differential equations in r, which suitable boundary conditions at the centre and 
surface, which can be solved numerically, with the frequencies LOnim as eigenvalues. How- 
ever, a great deal of insight into the properties of the solution and the relation of the 
frequencies to the structure of t h e sta r can be obtained from asymptotic analysis of 
the equations. iDeubner fc Goughl (|l984l ) derived an asymptotic equation, applicable for 
modes of high radial order, in terms of the quantity X — c^p^/^div5r, where Sr is the 
displacement vector: 



d^ 



^K{r)X 



where 



K 



1 



(2.7) 



(2.8) 



Thus K is controlled by three characteristic frequencies of the star: the Lamb frequency 
Si, with 



1(1 + iy 



(2.9) 
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the buoyancy frequency (or Brunt- Vdisdld frequency) N, 



Fi dr dr 

where g is the local gravitational acceleration, and the acoustic cut-off frequency oJc, 

where H = —{d\n p/dr)^^ is the density scale height. 

Equation (|2.7p determines the local properties of the eigenfunction. In regions where 
K{r) > the solution locally oscillates as a function of r, while the behaviour is exponen- 
tially increasing or decreasing where K{r) < 0. The intermediate points, where K{r) — 0, 
are called turning points. Often there is an interval [ri,r2] where K(r) > 0, with K be- 
ing negative just outside and the solution decreasing exponentially in the direction away 
from that interval. In that case the mode is said to be trapped between ri and r2 ; from 



JWKB analysis fe.g.. iGoughl 120071 ) one finds that the frequency approximately satisfies 
the dispersion relation 

K^/^dr ^ {n~l/2)T: . (2.12) 



ri 



The behaviour of the Lamb frequency and the acoustic cut-off frequency is relatively 
simple. The Lamb frequency generally decreases with increasing r||| for the low degrees 
that are relevant to stellar observations it is small near the surface. Here, typically, 
dominates in Eq. (12. 8p . In the atmosphere, where the temperature is approximately 
constant, H can be approximated by the (constant) pressure scale height and we obtain 

OJ, ~ L^catm = IgV^lp/P « MR-'^T;^'^ , (2.13) 

evaluated in the atmosphere; in the last proportionality we assumed the ideal gas law (see 
Eq. 12. 3p and that the atmospheric temperature is proportional to the effective tempera- 
ture Tcff . Below the surface Wc decreases rapidly with increasing depth. When uj < Wc.atm 
the solution decreases with height in the atmosphere; this corresponds to reflection of the 
mode at the surface so that it is trapped in the stellar interior. Waves with cj > Wc,atm can 
propagate out through the atmosphere and hence lose energy, leading to strong damp- 
ing. Thus Wc,atm defines the upper limit in frequency to trapped modes, at least in the 
adiabatic approximation. 

The behaviour of the buoyancy frequency is more complex. This is seen most clearly 
by assuming the ideal gas approximation, Eq. (j2.3p . Then we obtain, from Eq. (|2.10p . 
that 

=ff'-(Vad- V + V,J , (2.14) 
P 

where, following the usual convention (see Kawaler, this volume), 

^ _ dlnT ^ _ f dlnT \ ^ _ dln^ 

dlnp ' ^'^ \91np / ' ^ dlnp 

Since typically p increases towards the centre and hence with increasing p, as a result of 
nuclear burning, the term in typically gives a positive contribution to N"^ in the deep 
interior of the star. It should also be noted that the condition for convective instabilitjij 

f possibly with the exception of discontinuities in composition, caused by nuclear burning in 
the stellar core, which leads to discontinuities in density and hence sound speed 

J strictly speaking the Ledoux criterion 
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Figure 1. Evolution track for a IMq model, approximately calibrated to the Sun at the 

present solar age. 

is that N'^ < 0. In convection zones the composition is uniform and ~ 0; also, except 
near the surface V is nearly adiabatic and N'^ is only slightly negative. 

Solar-like oscillations are believed to be excited stochastically by the vigorous motion in 
the near-surface layers of convective envelopes (see Section [3] below). Thus stars showing 
such oscillations have outer convection zones and hence the buoyancy frequency is positive 
only in the deeper interior. Here we illustrate the characteristic frequencies and their 
dependence on the evolutionary state by considering selected models in a 1 Mq evolution 
sequence, illustrated in an HR diagram in Fig. [TJ A typical example of a rather unevolved 
star, at roughly solar age, is shown in Fig. [5] Here there is a fairly clear separation between 
Si and N. As a result there are essentially two different scenarios for having a positive K, 
also illustrated in the figure: by having uj > Si (and hence typically w ^ N) or by having 
oj < N (and hence typically oj <^ Si). The former case corresponds to acoustic modes 
(or p modes), where the restoring force is predominantly pressure, while the latter case 
corresponds to internal gravity waves (or g modes), where the restoring force is buoyancy. 

In discussing the properties of these modes in more detail, we neglect Uc in Eq. (|2.8p . 
except insofar as to use that it causes the required reflection at the surface of the p 



6 



J. Christensen-Dalsgaard: Asteroseismology of red giants 




Figure 2. Characteristic cyclic frequencies for Model 1 in the IMq evolution sequence illustrated 
in Fig.[T] The solid curve shows N/2tt, where A'' is the buoyancy frequency (cf. Eq. I2.10|) . whereas 
the dashed, dot-dashed and triple-dot-dashed curves show the Lamb frequency Si/2-n (cf. Eq. I2.9|I 
for I — 1, 2 and 3, respectively. The horizontal solid lines show typical propagation regions for 
a p mode (with / = 1, at high frequency) and a g mode (at low frequency). 



modes. For the p modes, assuming that ui ^ N, we then obtain from Eq. ()2.8|) that 

(2-16) 

with a lower turning point rt defined by S*; = a; or 

(2.17) 



c(rt) 



L 



where =1(1 
ranged to yield 



1). Also, the corresponding approximation to Eq. (|2.12l) can be rear- 

{n + a)Tr 



R 



1 



1/2 



dr 
c 



(2.18) 



Here we replaced the —1/2 in Eq. (I2.12p by a, which generally is a function of frequency, 
to account for the phase change in the reflection at the surface. Equation (|2.18l) has 
played a major role in the an a lysis o f solar oscillation frequencies; a relation of this form 
was first obtained by iDuvalll (Il982l) from analysis of observed frequencies. It explicitly 
demonstrates that in this limit the frequencies are determined by the sound speed and 
by a. 

In unevolved stars rt is close to the centre where c varies relatively slowly. Equa- 
tion (I2.18|) can then be expanded (|Goughlll986l . Il993f) to yield, in terms of the cyclic 
frequencies Vni = w„i/27r. 



- dn 



(2.19) 
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where 

A.j2p-r.)\ (2.20) 



and we replaced L by Z + 1/2 (see also Vandakurov 1967 : Tassoull 198d . ll990[ ) . Neglecting 



the small correction term c?„/ this shows that the frequencies are uniformly spaced in 
radial order n, with a spacing given by the large separation Av and such that there is de- 
generacy between Vni and Vn-i i+2- This degeneracy is lifted by dni which asymptotically, 
for main-sequence models, can be approximated by 

.„,.--4^.»r^£i. (2.21) 

47r^t'„; Jq dr r 

To this approximation (i„/, and hence the small separations Siyii^2{n-) = Vni — 
are very sensitive to the sound-speed gradient in the stellar core and hence, according 
to Eq. (|2.4I) . to the composition profile that has resulted from the nuclear burning. 
Thus for main-sequence stars the small frequency separations provide a measure of the 
evolutionary stage and hence the age of the star. Interestingly, although the derivation 
sketched above and the expression for dni do not hold for highly evolved stars such as 
red giants, the frequency pattern reflected in Eq. (|2.19p is still found (see also Section IH]). 
In the opposite extreme, where lo <^ Si and lu < N, we approximate K by 

With a relatively simple behaviour of N such as illustrated in Fig. [5] there are typically 
just two turning points ri, r2 where io — N, and Eq. (12.121) yields 



r2 



ri 



\ dr 



UJ^ J r 



L — = (n-l/2)7r. (2.23) 



li to <^ N almost everywhere on [ri,r2] we can approximate Eq. (|2.23p further by ne- 
glecting 1 compared with N'^/ui^, to obtain finally an approximate expression for the 
period 11 — 2ii juj: 

n= ^^5^(n-|-ag) , (2.24) 
Li 



where 

Ho = 27r2 ( / N— \ , (2.25) 



'■^ dr 



r 



assuming that iV^ > on [ri,r2], with = at ri and r2 (e.g., Tassoull 19801 ) . Also, 
ag, which may depend on refiects the actual phase shift at the turning points. Equa- 
tion (|2.24p shows that in this case we get periods that are uniformly spaced in the radial 
order n, with a spacing that is inversely proportional to L. 

The dependence of the characteristic frequencies and hence the oscillation properties 
on stellar parameters can to a large extent be characterized by simple scaling relations. 
An important example is the scaling of the acoustic cut-off frequency in Eq. (j2.13l) . Under 
many circumstances the structure of a star, as a function of the fractional radius x = r/R, 
appro ximately satisfies so-called homology scaling relations fe.g.. lKippenhahn fc WeigertI 
1993)0 

According to these pressure scales as GM^/R^, G being the gravitational con- 



stant, and density scales as M/R^. It follows that the squared sound speed c scales 



f A simple example, where the scalings are exact, are polytropic models. 
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Figure 3. Dimensionless characteristic frequencies (in units of {GM / B?Y^^) for models along 
the IMq evolutionary track illustrated in Fig. [T] Solid, dashed and dot-dashed curves are for 
Models 1, 2 and 3, respectively. The thick curves show the buoyancy frequency A'^ while the 
thinner curves show the Lamb frequency Si, for I = 2. The horizontal lines at the right edge of 
the plot similarly show the dimensionless value of the acoustic cut-off frequency (cf. Eq. I2.13P 
in the atmosphere of the models. 



as GM/R and the squared characteristic frequencies scale as GM/E?. This motivates 
introducing the dimensionless scaled frequencies by 



S 



GM 



o2 



GM' 



Col 



GM^ 



as well as (? = {R/ GM)c? . Then the asymptotic equations (|2.7p and 
pressed as 



da;^ 



-K{r)X 



where 



K 



in terms of the dimensionless frequency 



GM^ 



(2.26) 
can be ex- 

(2.27) 
(2.28) 

(2.29) 



With similar scalings the full equations of adiabatic oscillations can be expressed in a 
form that is homologously invariant. It follows that homologous models have the same 
dimensionless frequencies tT„i, and that the actual frequencies satisfy the scaling 



^„iOc(A//i?3) 1/2 0,(^)1/2^ 



(2.30) 



where (p) is the mean density of the star. The same scaling obviously applies, e.g., to the 
large frequency separation Aj/. 

Actual stellar models are never strictly homologous, and in particular the changing 
structure of the core as nuclear burning proceeds gives strong departures from these 
scalings. On the other hand, they are approximately satisfied in the outer layers which 
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predominantly determine the frequencies of acoustic modes, and hence for such modes 
the scahng in Eq. p.30p is a reasonable approximation. To illustrate this we consider the 
evolution of the frequencies along a 1 Af© evolution sequence, illustrated in Fig. [T] The 
dimensionless frequencies for Models 1, 2 and 3 are shown in Fig. [31 Evidently there is 
relatively little change in 52- On the other hand, the changes in N are dramatic, showing 
very strong departures from the homologous scaling. The increase in the depth of the 
outer convection zone, where \N\ ~ 0, is evident. Also, N increases strongly in the core. 
This is mainly the result of the increase in the gravitational acceleration (cf. Eq. 12.141) 
as the star develops a very compact helium core. Thus Model 3 has a helium core with 
a mass of 0.18 M©, contained within the inner 0.7 % of the stellar radius. In this model 
the convective envelope is so deep that it extends into a region where the composition 
has been changed by hydrogen fusion. This causes a discontinuity in the composition 
and hence the density, giving rise to a spike in the buoyancy frequency, as shown in the 
figure. 

In Fig. [3] are also indicated the dimensionless values Wc.atm of the acoustic cut-off 
frequency in the atmospheres of the three models. Since this is determined by the atmo- 
spheric properties of the star it does not follow the homologous frequency scaling. In fact, 
using Eq. (I2.13|) for i^c,atm — Wc,atm/27r and the scaling in Eq. (|2.30|) for Ai/ it follows 
that 



Since fc,atm is an upper limit to the frequencies of modes trapped in the stellar interior, 
this ratio, according to Eq. (12.191) . provides a measure of the number of acoustic modes 
that are trapped. This obviously decreases as the star evolves up the red-giant branch. 

Figure [3] also shows that in Models 2 and 3 the maximum of N in the core is substan- 
tially higher than Wc,atm- This means that all trapped nonradial modes satisfy lu < Si 
in the core and hence behave like g modes there, while they may behave as acoustic 
modes in the outer parts of the star where u! > N,Si. The overall character of the 
mode then depends on whether it has the largest amplitude in the g-mode or the p- 
mode region, and hence, effectively, on the behaviour in the intermediate region. If the 
eigenfunction decreases exponentially with depth in this region the mode has its largest 
amplitude in the outer region of the star and it may be said to be p- dominated. In the 
opposite case, with an eigenfunction increasing with depth, the amplitude is largest in 
the core and the mode is g-dominated. These two conditions essentially require the fre- 
quency to be determined such that the mode resonates with the outer acoustic, and 
the inner buoyancy-driven, cavity, respectively. For the p-dominated modes this leads to 
Eq. (|2.18|) for the frequencies, which may be cast, at least approximately, in the form 
of Eq. (j2.19|) . For the g-dominated modes we similarly find that the frequencies satisfy 
Eqs (|2?23| or (|2?24l) . For evolved red giants like Model 3, with a very high buoyancy fre- 
quency in the core, the period spacing Hq/L (cf. Eq l2.25p becomes very small and hence 
the density of the g-dominated modes is much higher than the density of p-dominated 
modes, whose frequency spacing is given by Av. It should also be noted that the high 
density of g-dominated modes correspond to an extremely high number of nodes in the 
eigenfunctions in the core, requiring substantial care in the numerical computation of the 
oscillations. 

The properties of the modes can be illustrated by considering the normalized mode 
inertia 



^c,atm 



(2.31) 



Ai^ 



E = 



IvP\Sr\'dV 
M\Sr\l 



(2.32) 
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Figure 4. Dimensionless inertia (cf. Eq I2.32p for Model 3 in the 1 Mq evolution sequence 
illustrated in Fig. [T] Radial modes are shown with circles connected by solid lines, / = 1 modes 
with triangles connected by dashed lines and I — 2 modes with squares connected by dot-dashed 
lines. The vertical dotted lines mark the / — 2 modes for which the eigenf unctions are illustrated 
in Fig. O 

where the integral is over the volume V of the star and |5r|g is the surface value of 
the squared norm of the displacement vector Sr. This is shown in Fig. 2] for modes of 
degree / = — 2 in Model 3. The radial modes {I = 0) are purely acoustic, with an 
inertia that decreases with increasing frequency. For I = 1 and 2, on the other hand, 
there is a dense set of g-dominated modes, with acoustic resonances where the inertia 
approaches the value for the radial modes at the corresponding frequency. The pattern of 
these resonances is very similar to Eq. (|2.19p , with the I = 1 resonances roughly half-way 
between the I = frequencies and a small separation between an ^ = 2 resonance and 
the neighbouring / = mode. 

Figure [5] shows a few selected eigenfunctions, in the form of scaled radial component 
of the displacement, normalized to unity at the surface. The comparison between the 
neighbouring I = and 2 modes in panel a) shows that they are very similar in the 
bulk of the star where the I — 2 mode is predominantly of acoustic character. In the 
core the eigenfunction for I = 2 varies rapidly but at generally small amplitude for 
this p-doniinated mode, as seen in panel b). The behaviour in the core for the mode 
shown in panel c) is very similar, but at a much higher amplitude relative to the surface 
normalization, and corresponding to the local maximum in the inertia. We discuss the 
consequences for the observable visibility of the modes in Section [5] below. 

3. Damping and excitation of solar-like oscillations 

When nonadiabatic effects are taken into account, the oscillation frequencies are no 
longer real. Assuming again a time dependence as exp(— iwi) and writing w = -f iwj in 
terms of real and imaginary parts, the real part Wr determines the frequency of the mode 
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Figure 5. Scaled eigenfunctions, normalized at the surface, for Model 3 in the 1 Mq evolution 
sequence shown in Fig.[T] The quantity shown is rp^^^^r, where is the amplitude of the vertical 
displacement. In panels a) and b) the dashed line shows a radial mode with cyclic frequency 
v = 219.6 /iHz and mode inertia (cf. Eq. I2.32[) E = 2.35 x 10^^, and the solid line shows a 
p-dominated / = 2 mode with u = 218.1 /^Hz and E — 3.12 x 10~^. In panel c) is illustrated a 
neighbouring g-dominated I = 2 mode with u = 226.1 /iHz and E = 1.00 x 10~*. 



whereas the imaginary part uJi determines the growth or decay of the amplitude, Wi < 
corresponding to a damped mode. The determination of the damping or excitation of the 
modes requires a solution of the full nonadiabatic oscillation equations, introducing also 
the energy equation and equations for the energy generation and transport. The latter 
is given by the perturbations to the flux F = Fj-ad + Fcon which, as indicated, has both 
a radiative (Frad) and convective (Fcon) contribution. The radiative part can be dealt 
with in a relatively straightforward manner, while the treatment of the convective part 
remains highly uncertain. Furthermore, a significant contribution to the damping may 
come from the dynamical effects of convection through the perturbation to the Reynolds 
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stresses, often approximated by a turbulent pressure pt, which again is uncertain (e.g., 
lHoudek.2010al 

A convenient way to discuss the damping or excitation of the modes is to express the 
growth rate, approximately in terms of the work integral W = Wg + Wt, as 



J 



J 



(3.1) 



with 





= Re 


[X 


p 


Wt 


= -Wi 


Im 


[X 


J 


= 2^2 


Jv 



-(Fa - l)<5(pe-divF)dT/ 
P 

^dV = 2uj'iM\5rJ^E . 



(3.2) 



Here 5 denotes the Lagrangian perturbation, i.e., the perturbation following the motion, 
Fa — 1 = ((?lnT/i91np)ad, and e is the rate of energy generation per unit mass; also the 
star indicates the complex conjugate. 

Equations (j3.ip and p.2p have a simple physical meaning: In the expression for VVg 
6{pe — div F) gives the perturbation to the heating rate; thus gets a positive con- 
tribution, and hence a contribution to the excitation of the mode, where heating and 
compression, described by 5p/ p, are in phase. Heating at compression is the basis for 
the operation of any heat engine, and hence modes that are unstable because of this can 
be said to be driven by the heat-engine mechanism^ Depending on the phase relation 
between compression and 5pt the term in Wt may contribute to the damping or driving 
of the mode. 

As already indicated, the treatment of the convection-pulsation interaction is a ma- 
jor uncertainty in nonadiabatic calculations for stars with significant outer convection 
zones, where the convective timescales and the pulsation periods are typically similar. 
Two generalizations of mixing- length theory have seen fairly widespread use. In one, 
based on an original model by lUnno ( 1967 ). steady convective eddies are considered, 
wi th a bal a nce b etween buoyancy and turbulent viscous drag. In the second, developed 
bv iGougii (1977), convective eddies are continually created and destroyed. The effects 
of the perturbations associated with the pulsations on these physical models can then 
be described, leading to expressions for the perturbations to the co nvective flux and 
Reyn olds stresses. Unno's description has been further developed (see iGrigahcene et al.l 



2005, and references t herei n), while Gough's de scription was developed and applied by. 



BalmforthI (1992) and Houdek et al. ( 19991 ). With judicious (but not unreasonable) 



e.g. 

choices of parameters both formulations can predict the transition to stability at the red 
edge of the Cepheid instability strip, as well as yield results for solar models that are 
consistent with the observationally inferred damping ratesU Interestingly, the dominant 
mechanism resulting in damping of the modes differs between the two formulations: in 
the calculations based on Gough's formulation damping is dominated by the perturbation 
to the turbulent pressure, while the results based on Unno's formulation show damping 



f Since the perturbation to the heating is often determined by the perturbation to the opacity 
K this is also known as the k me chanism. 

X A different formulation by IXiong et al.l (|l997l ). based on equations for the second- and 
third-order correlations, similarly predicts the red edge but apparently finds that some of the 
observed solar modes are unstable (jXiong fc Dengif20ldl ). 
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Figure 6. Power spectrum of a solar radial mode, from 8 years of BiSON observations. The 
white curve shows a Lorentzian fit (cf. Eq. 13. 3p to the observed spectrum, while the smooth grey 
curve shows the sam e fit, but multiplied by a factor three. Data courtesy of W. J. Chaplin; see 
IChaplin et all (|2002h . 



through the perturbation to the convective flux. A more detailed comparison of these 
two sets of results is clearly called for. 

The cause of the solar oscillations is now almost universally thought to be stochastic 
excitation by near-surface convection of the intrinsically stabl e modes. Sto chastic exci- 
tation of waves in the solar atmosphere was first coi isidered by SteinI ( 19681), fo llowed by 
the application to the excitation of normal modes bv lGoldreich fc KeelevI (fl977h. The ba- 
sic pro perties of damped oscillations excited by tur bulence was considered bv iBatcheloj 
( 1956[ ) (see also Christensen-Dalsgaard et al. 1989t ). The result is that the power spec- 
trum of a single mode, of natural frequency and damping rate Wi, is 



Pic) 



(3.3) 



where Pi{uj) is the power spectrum of turbulent forcing. Assuming that Pf varies slowly 
with frequency the result is a stochastic function modulated by a Lorentzian envelope, 
with a full width at half maximum of 2|aji|. An example, from BiSON observations 
of a solar radial mode, is shown in Fig. [51 From the fitted Lorentzian it is obviously 
possible to determine the damping rate jwij. Further analyses of the statistica l prop- 
erti es of the stochast i cally excited oscillators were carried out by Kumar et al.1 ( 19881) 
and IChang fc Gough (Eiii). It was shown by IChaplin et aP (11997^ that the observed 
distribution of solar oscillation amplitudes was consistent with the predictions of these 
analyses. 

An e arly rough prediction of the excitation of so lar-like modes in o ther stars was carried 



out by [Christensen-Dalsgaard fc Frandsen ( 19831 ). It was shown bv iKieldsen fc Bedding 
(119951 ) that the resulting velocity amplitudes approximately scaled as L/M, where L is the 



surface luminosity o f the star. Later, more careful investigations (e.g.. lHoudek et alll 19991 : 
Samadi et al.l[2007l ) confirm that the amplitudes increase with L/M although generally 
to a power slightly less than one. 

The excitation of a specific mode clearly depends on the detailed property of the mode, 
the amplitude being determined by the balance between the stochastic energy input from 
convection and the damping. The energy input increases as a high power of the convective 
velocity (Stein. 1968<) and hence is dominated by the near-surface layers where convection 
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is most vigorous. Here the properties of the oscillations depend predominantly on the 
frequency, with little dependence on the degree (at least for the low-degree modes that 
are relevant in distant stars). The o utcome is that the mean square amplitude can be 
written as fe.g.. IChaplin et al.ll2005l) 



1 Vfjuj) 
E\oJi\ E 



(3.4) 



where A might be the surface velocity or the relative intensity fluctuation; here Vi{uj), 
which depends on freque ncy but in gene ral not on degree, describes the energy input 
from convection (see, e.g., HoudeklboiObl for a review). Also, the damping rate is given 
by Eqs (|3.ip and p.2|) . In the common case where the integrals in Wg and Wt are dom- 
inated by the near-surface layers it follows that \u!i\E is independent of E at constant 
frequency; consequently, in this case {A^) oc E~^ at constant frequency. This is impor- 
tant in the discussion of the excitation of solar-like oscillations in red giants where we see 
large variations in E over a narrow frequency range (cf. Fig. the g-dominated modes, 
with large inertias, are clearly expected to have a much smaller root-mean-square ampli- 
tude than the p-dominated modes with low inertia. However, the question of the mode 
visibility, and hence the likelihood to detect them, is more complicated; we return to this 
in Section m 

The balance between energy input and damping leads to a characteristic bell-shaped 
amplitude distribution of the mean amplitude, although with fluctuations around the 
mean reflecting the stochasti c nature of the excitation. The shape of the amplitude 
distribution was discussed by iGoldreich et al. ([l994) . They found that the increase in 



amplitude with frequency at low frequency was predominantly caused by the change 
in the shape of the eigenfunction as the upper reflection point, where lo = lo^ moves 
closer to the surface with increasing lo. At high frequency this upper reflection point is 
very close to the photosphere, and the amplitude decreases with increasing frequency 
due to an increasing mismatch between the oscillation periods and the timescales of 
the convective eddies that dominate the excitation. As the frequency approaches Wc,atm 
beginning energy loss through the stellar atmosphere may also contribute to the damping 
and hence further decrease the amplitudes. Indeed, the variation with frequency of the 
damping rate clearly also plays a role for the the variation in amplitude. 

It has been found observationally that the cyclic frequency t'max o f maximum ampli 



Brown et al. 1991 



tudc appears to scale with the aco ustic cut-off frequency Vc.atm (e.g., 
iKjcldsen & Bedding 1995; Bedding "fcKieldsenl [2003t IStello et al.ll2008t ) (see also Bed- 



ding, this volume), and hence follows the scaling in Eq. (|2.13l) . This scaling has been 
a very useful too l in characterizing th e observed properties of stars showing solar-like 
oscillations (e.g., iKallinger et al. 201Clt ). bu t the relation b e tween fmax and I'catm still 
defies a definite theoretical understanding. iBelkacem et al.l (|2011l ) noted that the loca- 
tion of J^max is closely related to the plateau in damping rates at intermediate frequency, 
apparently caused by a resonance bet ween the pulsa tion period and the local thermal 
timesc ale, a feature already noticed by iGougb ( 1980l ) and Christensen-Dalsgaard et al 
( 1989( ). However, this does not fully explain the observed tight relation. 



4. Do red giants have nonradial oscillations? 

As pointed out bv lKieldsen fc Beddind ()l995l) the amplitude predictions bv lChristensen-Dalsgaard fc FrandsenI 
(|l983l ) scaled as L/M. Thus one could expect to see solar-like oscillations in red giants 
with q uite substantial amplitudes. Early e videnc e was found for such oscillatio ns in Arc- 
turus (jSmith et al.lll987l: llnnis et al.lll988l) while lEdmonds fc Gillila^ (|l996h identified 
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Figure 7. Power spectrum of the G7 giant ^ H ydrae, from one month of Doppler velocity 
measurements. Adapted from iFrandsen et aLl (|2002l ). 



oscillations in K giants in the globular cluster 47 Tuc. lChristensen-Dalsgaard et al.l (|200lh 
noted that the statistics of very long-period semi-regular variables observed by the Amer- 
ical Association of Variable Star Observers appeared to obey amplitude statistics con- 
sistent with stochastic excitation, strongly suggesting that the variability of these stars 
was also solar-like, although at month-long periods and with amplitudes allowing simple 
visual observations. Photometric observations of the oscillations in Arcturus were carried 
out by Retter et al.l (|2003) with the WIRE sat ellite. Analyses of lar ge samples of red gi- 
ants observed by OGLeQ fe.g. lKiss fc Beddinel[200l |2004 ISoszvhski et al.ll2Q07h showed 
oscillations in several modes, and Dziembowski fc SoszyhskT i 20101 ) confirmed that the 
most likely cause was solar-like oscillations. In a further analysis of ground- based survey 



observ ations of red giants in the Galaxy and the Large Magellanic Could iTabur et al 



( 2010h demonstrated a continuous transition between pulsations in G and K giants and 
those seen in very luminous M giants. 

The first detailed o bservations of the oscillations of a red giant was carried out by 
Frandsen etHI (|2002[ ) for the G7 giant ^ Hya. The resulting power spectrum is shown 
in Fig. [71 This clearly shows peaks at uniformly spaced frequency, in accordance with 
Eq. (|2.19p . and a power envelope very similar to what is observed in the Sun, strongly 
supporting the solar-like nature of the oscillations. From analysis of the spectrum the 
separation between the peaks was determined to be ~ 7 ^Hz. 

Modelling of ^ Hya to match the well determined luminosity and effective temperature 
is shown in Fig. [H] Stars ascend the red-giant branch while burning hydrogen in a shell 
around an inert and almost isothermal helium core. As the mass of the core increases, so 
does the luminosity and the temperature of the core. At the tip of the red-giant branch 
the temperature reaches a sufficiently high value, around 10® K, where helium fusion 
becomes significant. The resulting change in structure, with an expansion of the core, 
leads to a reduction in the surface radius and luminosity, with the star entering a phase 
of stable helium burning. Even during this phase, however, a substantial part of the total 
energy production takes place in the hydrogen-burning shell. As the phase is relatively 
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Figure 8. Evolution track for a model of ^ Hydrae, with mass M = 2.8i\fQ and heavy-element 
abundance Z = 0.019. The plusses are placed at 5 Myr intervals along the track. The box marks 
the 1—a error box for the observed effective temperature and luminosity. 



extended, with a duration of around 20 per cent of that of the central hydrogen-burning 
main sequence, the chances of observing a star in this phase are fairly high; in stellar 
clusters such s tars populate the so-c alled 'red clump'. 

As noted bv lTeixeira et al. I (|2003l ) C Hya could be identified as being in the ascending 
phase on the red giant branch. However, since this phase is extremely brief a more likely 
identification would be with the core helium burning phase. In either case, a match to the 
observed frequencies was possible only if the observed frequency spacing was identified 
with Al/, such that only modes of one degree, presumably radial modes, were observed. 
In contrast, from a strict application of Eq. (|2.19p one would have expected the spacing 
to be Ai//2, and hence Av ~ 14/iHz. However, given the scaling in Eq. (|2.30p this would 
be entirely inconsistent with the loc ation of the star in the HR di agram. 

To account for this identification IChristensen-Dalsgaar 3 (|2004l) considered the possi- 
bility that all nonradial modes in red giants would likely be heavily damped and hence 
have small amplitudes. He noted that Wg (cf. Eq. 13. 2p depends on (5(divF), i.e., on the 
derivative of SF. In turn, roughly SF cx V((5r) and 6T/T oc 5p/p oc div Sr. It follows 
that Wg depends on the third derivative of Sr. Thus even if the amplitude of Sr in the 
core is quite small, as is the case for p-dominated modes (cf. Fig. [5]), (5(divF) could be 
quite large for a mode varying very rapidly as a function of r, leading to a substantial 
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Figure 9. Power spectra of solar-like oscillations in red giants, from 5 months of observations 
with the CoRoT satellite. The panels are labelled by the CoRoT identification number. Radius 
and lu minosity increase towards the top. Figure courtesy of J. De Ridder; from lDe Ridder et al.l 
l|2009l ). 



contribution to the damping from the core for all nonradi al modes. Although the anal- 
ysis was entirely qualitative. IChristensen-Dalseaardl ( 2004) found some support for this 
conjecture in the analyses by Dziembowskil ( 19771 ) and Dziembowski et al. ( 2001 ). 

This pessimistic expectation for red-giant seismology has been complete l y over turned 
by observations and an improved theoretical understanding. Hekker et al. ( 2006[ ) found 
evidence for nonradial pulsation in the line-prof ile variation of a f ew re d giants showing 
solar-like oscillations. Even more dramatically, iDe Ridder et al. I (l2009l) clearly demon- 
strated the presence of nonradial solar-like oscillations in a large number of red giants 
from early CoRoT observations. Figure [5] shows resulting power spectra. The power en- 
hancement from the solar-like oscillations is evident; as stars on the red-giant branch have 
roughly the same effective temperature, the frequency at maximum power decreases with 
increasing radius and luminosity in accordance with Eq. (j2.13p . The observed frequencies 
satisfy a relation similar to Eq. (I2.19p . with definite evidence for nonradial modes. This 
is most clearly illustrated in an echelle diagram (see Bedding, this volume), where the 
frequencies are shown against the frequency reduced modulo Av, effectively dividing the 
spectrum into segments of length Av and stacking the segments. This has been done for 
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Figure 10. Echelle diagram of the star 101034881 in Fig. |9] This corresponds to stacking shces 
of the power spectrum of a length given the large separation Av — 4.81 /iHz: the abscissa is the 
frequ ency modulo Av whereas the ordinate is the frequency. Figure courtesy of J. De Ridder; 
from lDe Ridder et al.1 (|2009l ). 



one of the stars in Fig. [TUl In accordance with Eq. (|2.19|) the three columns of points 
correspond to modes of degree I = and 2 to the right, separated by the small separation, 
and modes with I = 1 to the left. Analyses o f red giants observed by CoRoT and Kepler 
(e.g., Huber et al. 20101 : Mosser et al. 2011a ) have shown that this frequency pattern is 
universal amongst red giants. 

The very long time series obtained by CoRoT and in particular Kepler have allowed 
fits to the power spectra in terms of the Lorentzian profiles (cf. Eq. 13.31) and hence a 
determination of the mode damping rates or, equivalently, the lifetimes t = l /lcJil. For red 
giants this yields lifetime s near the maximum amplitude of around 15 days ( Huber et al.l 
120101: Baudin et al. '2011V Interestingly, this is very similar to the value obtained for ^ 
byl 



Hya bv lHoudek &: Gough (2002), using a convection treatment based on iGough (IQTtI ). 



5. Amplitudes of solar-like oscillations in red giants 

To interpret the observations a more detailed analysis of the mode excitation and visi- 
bility is required. The analysis leading to Eq. p.4|) determines the mean square amplitude 
of the mode, corresponding essentially to t he area under the co rresponding peak in the 
power spectrum. However, it was noted by Chaplin et al.l ( 2005f) that ability to detect a 



mode in the power spectrum depends instead on the peak height H, related to by 

= \uji\H , (5.1) 

smce \uii\ determines the width of the peak. If the damping is dominated by the near- 
surface layers, it was argued above that both and \oJi\ are proportional to E~-^ at 
a given frequency. It then follows from Eq. (|5.ip that H is independent of E, i.e., that 
all modes in a given frequency interval are excited to roughly the same amplitude. For 
the dense spectrum of modes in a red giant, such as illustrated in Fig. [H the result 
would be a very dense observed power spectrum, with no indication of the acoustic 
resonances approxir nately satisfying Eg. (I2.19|) . This is in obvious contradiction to, e.g., 
the observations bv lDe Ridder et all (|2009i) . 



This argument, however, neglects the fact that the oscillations are observed for a finite 



f for simplicity, we drop (...) in the following 
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Figure 1 1 . Evolution tracks for models of mass 2 and 3 Mq , for which iDupret et al.l l|2009D 
considered the stochastic excitation of solar-like oscillations. Figure courtesy of M.-A. Dupret. 



time T which, for the g-dominated modes, is typically shorter than the lifetime r of the 
mode. The resulting peaks in the power-density spectrum has a width that, for T ^ t, 
is proportional to and hence a correspondingly smaller peak height, very likely 

making the modes i nvisible. As a convenient interpolation between the two extremes. 



Fletcher et al.l (|2006[ ) proposed the relation 



42 

H(x— — — — . (5.2) 

l^il + 2/r ^ ' 

A detailed and very illum i nating analysis of the excitation of modes in red giants was 
carried out bv lDupret et al. They considered several red-giant models, with loca- 

tions in the HR diagram as illustrated in Fig. 1111 carrying out nonadiabatic calculations 
of nonradial modes using the convection for malism of Grigahccnc ct al. (2005). From a 



model of the stochastic energy input (e.g., Belkacem et all feoOGl) they then computed 



the expected amplitudes of the modes and the resulting power-density diagrams. They 
assumed an observing period T of 150 days, corresponding to CoRoT long runs. Rather 
than using Eq. (lOt they effectively adopted H cx tA^ for r > T /2 and iJ cx T 
otherwise, with a continuous transition between the two cases. Here we briefly consider 
their results for Models A and B. 

Figure [12] shows the mode lifetimes for modes of degree Z = — 2 in Model A. In this 
case the damping is dominated by the near-surface layers and hence, as is obvious in the 
figure, the variation of r reflects the behaviour of E which is qualitatively similar to the 
case illustrated in Fig. HI For I — 2 the g-dominated modes have lifetimes substantially 
higher than the assumed T/2, indicated in the figure by a horizontal dashed line. On 
the other hand, the I = 1 modes at higher frequency all have lifetimes at or below the 
observing time, and hence would be expected to reach similar peak heights in a given 
frequency interval. This is confirmed by Fig. [T3] which shows the predicted peak heights. 
The overall envelope seems to be somewhat skewed towards lower frequency, relative to 
what is expected from the observed scaling with the acoustic cut-off frequency; given 
the uncertainties in modelling both the damping and the energy input, this is hardly 
surprising. However, it is striking that while the I = 2 modes show strong acoustic 
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Figure 12. Computed lifetimes of Model A in Fig. [11] at the base of the 2Mq red-giant branch. 
Modes of degree / = 0, 1 and 2 are shown, respectively, by plusses connected by solid lines, open 
circles connected by dashed lines and closed circles connected by dotted lines. The horizontal 
dashed line marks a damping time of 75 d, corresponding to the border between unresolved 
and resolved modes for a CoRoT long run of 150 days. Figure courtesy of M.-A. Dupret; from 
iDupret et aT] l|2009l l. 
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Figure 13. Estimated peak heights in the power spectrum of solar- like oscillations in Model A 
in Fig. 111! at the base of the red-giant branch; for symbol a nd line styles, see caption to Fig. 1121 
Figure courtesy of M.-A. Dupret; from lDupret et ahl (|2009f) . 



resonances, with small peak heights for the g-dominated modes, the / = 1 modes at high 
frequency all reach peak heights comparable with tho se for the radi al modes. It is obvious 
that such a spectrum can be quite complex (see also lDi Mauro et a l. 2011). 

The situation is rather different for the more evolved Model B. As illustrated in Fig. [Ml 
there is now strong damping in the core for the g-dominated modes with I = 2; for I = 1 
the same is the case at relatively low frequency, while the surface layers dominate at high 
frequency. The damping in the core is a result of the very rapidly oscillating eigenfunctions 
in the core, as discussed above. Owing to the increased central condensation the buoyancy 
frequency is substantially higher than in Model A, leading to a higher vertical wave 
number and hence stronger damping in the core. For p-dominated modes, on the other 
hand, the damping is still dominated by the near-surface layers. Together with a strong 
trapping of the g-dominated modes in the core, this leads to a much stronger contrast 
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Figure 14. Work integrals for the intermediate red-giant branch 2 A'/q Model B in Fig. [TT] this 
essentially corresponds to the partial integrals of Wg-j- Wt, defined in Eq. (I3.2|l . and here plotted 
against temperature T in the model. The top panel is for a p-dominated mode while the bottom 
pane l is for a g-dominate d mode largely trapped in the core. Figure courtesy of M.-A. Dupret; 
from lDupret et all l|2009D . 




Figure 15. Computed lifetimes for the intermediate red-giant branc h 2 Mp, Mode l B in 
Fig- [HI see caption to Fig. [12] Figure courtesy of M.-A. Dupret; from lDupret et al] (|2009l 'l. 



in lifetime between the p- and g-dominated modes, with the latter having lifetimes very 
substantially higher than the observing time (cf. Fig. [15]). The effect of this on the 
predicted power density spectrum is evident in Fig. 1161 there are now obvious acoustic 
resonances for both I = 1 and 2. The I = 2 resonances are confined to very narrow 
frequency intervals, caused also by the more extended evanescent region between the 
peak in the buoyancy frequency and the region of acoustic propagation, and typically 
involve only 1-2 modes at non- negligible peak heights, For I = 1, on the other hand. 
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Figure 16. Estimated peak heights in the power spectrum of solar-like oscillations in the inter- 
mediate red-giant branch 2 Mq Model B in Fig. 11 11 for sym bol and line styles, see caption to 
Fig. \n\ Figure courtesy of M.-A. Dupret; from iDupret et al.i pOOgP . 



there is typically a cluster of peaks with heights that may make the mo des observable. As 
disc ussed by Bedding (th is volume) this was found to be the case bv lBeck et al] (|201lh 
and iBedding et al. We discuss the significance of these results in the following 

section. 



6. Asteroseismic diagnostics in red giants 

Detection of solar-like oscillations in red giants, such as illustrated in Fig. [51 almost 
immediately provides a measure of the large frequency separation Av and the frequency 
Vmax of maximum power. These basic quantities already have a substantial diagnostic po- 
tential: Au approximately satisfies the scaling (Eq. I2.30p with the mean density and it can 
be assumed that fmax scales as the atmospherical acoustic cut-off frequency (Eq. I2.13P . 
If Toff is known from photometric or spectroscopic observations one can then determine 
the radius and mass of the star as 



R _ t^max / Av \ / Toff 
Rq 



1/2 



(6.1) 



and 



(6.2) 



where denotes solar values (e.g., Kallinger et al.l[201ol). Such analyses provide unique 



^ \ : 

possibilities for population studies of field red giants fe.g.. lMiglio et al.ll2009HMosser et al 

[2010; Hekkcr ct al. 2011). Stellar modelling provides further constraints on the relation 
between T^,ff, M and R and hence a more precise determination of stellar properties 



(|Gai et al.l 2011 ), as used in an a nalysis of red giants in two of the open clusters in the 
Kepler field bv lBasu et al.l (|201l[ ). 

As discussed in connection with Eq. (I2.2ip . for main-sequence stars the small frequency 
separations 5h'ii-f-2{n) = Vni — Vn-11+2 provide a measure of the age of the star. This no 
longer holds for red giants. Here in most cases 5vii+2 scales essentially as the large fre- 
quen cy separation and hence mainly reflects the variation of the m ean density of the star 
(e.g.- lBedding et al.llioiot iHuber et"aI1l2010HMosser et al.llioilal ). This is a natural con- 
sequence of the fact that the propagation region of the nonradial acoustic modes for such 
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Figure 17. The sohd curve shows the period spacing between adjacent I = 1 modes in Model 3 
in the IMq sequence in Fig. [1] The horizontal dashed line shows the asymptotic value for pure 
g modes, from Eq. (|6.4p . For comparison, the short-dashed curve shows the logarithm of the 
normalized inertia Qni (cf. Eq. 16. 3p on an arbitrary scale. 



stars lies in the convective envelope, which probably changes in a largely homologous fash- 
ion as the properties of the star change. A detailed analysis of the diagnostic potential of 
acoustic modes in red giants was carried out by Monta lban et a l. (2010). They also consid- 
ered other types of small frequency separations, such as 6voi{n) — (j^„o + t'n+i o)/2 — i^ni- 
Interestingly, they found that for models with somewhat shallower convective envelopes, 
particularly amongst core helium-burning stars, where the lower turning point of the I — 1 
modes was in the radiative region, the average {5vqi) / i^v showed a marked dependence 
on the distance from the base of the convective envelope to the turning point. 

The asymptotic analysis leading to Eq. (|2.19p for the acoustic-mode frequencies as- 
sumes that the equilibrium structure varies slowly compared with the eigenf unctions. 
There are regions in the star where this does not hold, with variations in the sound 
speed on a scale that is comparable with or smaller than the local wavelength. This 
causes periodic variations in the frequencies with properties that depend on the depth 
and nature of the sound-speed feature. In red giants the most important example of such 
an acoustic glitch is related to the second ionization of helium which causes a dip in Fi 
and hence (cf. Eq. 12. 2[) a localized variation in the so und speed. This variation provides a 
diagnosti cs of the helium abundance in the star (e.g..'Gough"l990':^Vor ontsov et al. 1991 



Basu et al . 2004; Mon tciro k, Thompson 2005; Houdck & Gough 2007). A clear signature 



of this effect was found by iMigho et all (I2OIOI1 in a red gian t observed by CoRoT. They 
noted, following an earlier discussion by iMazumdan (|2005l ) of the diagnostic potential 
of such glitches, that the inferred depth of the sound-speed feature together with large 
frequency separation provide a purely asteroseismic determination of the mass and radius 
of the star. This feature can probably be detected in a substantial fraction of the red 
giants observed with CoRoT and Kepler, promising a determination of the so far poorly 
known helium abundance in these stars and hence important constraints on the chemical 

evolution of the Galaxy. 

Prom the analysis bv iDupret et al. ( 20091) . discussed in Section [5l we noted that for 
I = 1 several modes might be excited to observable amplitudes in the vicinity of an acous- 
tic resonanc e. This is cons i stent with the rather larger scatter for these dipolar modes 
observed bv bedding et all (|2010l ) in an collapsed echelle diagram based on early Kepler 
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Figure 18. 'Hertzsprung-Russell diagram', in terms of the asymptotic dipolar g-mode spacing 
AIIi (cf. Eq. I6.4|l and luminosity, for the 1 Mq evolution track illustrated in Fig. [l] 



data and, as noted by iMontalban et"aL ( 2010() . has important diagnostic potential. As 
discus sed by Bedding (t his y olume) this pot e ntial was dramatically realized in the analy- 
ses by lBeck et al. ( 201ll) and Bedding et al. ( 2011 ). The groups of dipolar peaks near the 
acoustic resonances were found to haye a regular structure which allowed determination 
of well-defined period spacings, reminiscent of the behaviour of g modes (cf. Eq. I2.24p 
and matching a corresponding behaviour in the computed frequ encies of stellar mod els. 
Such behaviour has also been found in analysis of CoRoT data (^Moss er et al. I l2011bl) . 

To illustrate the properties of the computed period spacings Fig.lTTlshows the difference 
in period between adjacent dipolar modes in Model 3 in the 1 M© evolution sequence 
illustrated in Fig. [TJ It is clear that, particularly in the dense region of g-dominated 
modes at relatively low frequency, there is an almost uniform spacing, interspersed with 
decreases which are associated with the acoustic resonances. To illustrate this relation, 
the figure also shows log Qni where Qni is the mode inertia normalized by the inertia of 
a radial mode at the same frequency. Specifically, 



nt 



inl — 



Eo{uJnl) 



(6.3) 



where Eni is the inertia of the mode and Eq(uj„i) is the inertia of radial modes, inter- 
polated to the frequency uJni of the mode. It is obvious that the variation in the period 
spacing is closely related to the variation in the inertia. 
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Figure 19. Cyclic buoyancy frequencies N/2tv in Model 3 (solid), Model 4 (dashed) and 
Model 5 (dot-dashed) in the 1 Mq evolution sequence illustrated in Figs [1] and 1181 



The uniform period spacing is clearly only realized for the most g-domin ated modes 
whose high inertia makes detection unlikelyHI However, it was demonstrated bv lBedding et al 
(|201lh that with the Kepler data the pure g-mode period spacing can be determined from 
the observed frequencies in data of sufficient quality; this was illustrated in a 'g-mode 
echelle diagram', dividing the spectrum into segments of fixed length in period (see also 
Bedding, this volume). In other cases the period spacings can be inferred fr om analysis of 
the po w er spectra and e xtrap olated to the value for pure g modes (see also lMosser et al] 
2011bh . iBedding et al.l (j2011l ) showed that the inferred period spacing is substantially 
higher for clump stars in the core helium burning phase than for stars on the ascending 
red-giant branch, providing a clear separation into these two groups of stars that are 
superficially very similar. 

From the asymptotic results in Eqs (|2.24[) and (|2.25p the period spacing for pure dipolar 
g modes is approximately given by 



N 



dr 



(6.4) 



This is illustrated by the dashed horizontal line in Fig. [171 it is clearly in good agreement 
with the value obtained from the most strongly trapped modes. Given the AIIi, as 
discussed above, can be obtained from the observations, we therefore obtain diagnostics 
of the buoyancy frequency in the core of the star. 

This can be used to illustrate the evolution of the star in an 'HR' diagram, plotting 
L against AIIi, as done in Fig. [181 It is evident that AIIi varies little as the star moves 
up the red-giant branch. This behaviour can be understood in terms of the variation in 
the buoyancy frequency, illustrated for three representative models in Fig. 1191 On this 
logarithmic r scale the main effect is a shift towards smaller radii as the core contracts, 
with little change in the shape of N. Noting that the integral in Eq. (|6.4p is in terms of 
logr it follows that there is little change in AHi with evolution. It should also be noticed 
that A'' is depressed in the core of the star. This is probably caused by the increasing 
electron degeneracy of the gas, leading to dlnp/dlnp approaching 5/3 



I 



f except perhaps in extremely long timeseries, exceeding the lifetimes of even these modes 
J In white dwarfs this leads to a very small buoyancy frequency in the interior of the star, 
confining the g modes to the outer layers (|Fontaine fc Brassard. 2008i l. 
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Figure 20. Evolution track for a 2.5 Mq model. The track goes through helium ignition at the 
tip of the red-giant branch and ends near central helium exhaustion. 

As discussed by Bedding (this volume) there is a striking difference in AIIi between 
stars in the hydrogen shell-burning phase, as the 1 Mq model illustrated above, and 
models that burn helium in the core. Unfortunately ASTEC does not allow computation 
of a model through a helium flash, and hence the evolution of the 1 Mq model cannot be 
followed beyond helium ignition. For more massive stars the helium ignition takes place 
in a more quiet manner which can be followed by ASTEC. Figure [20] shows the evolution 
track of a 2.5 Mq model up to a point near exhaustion of helium at the centre. 

In the ( AIIi , L) diagram this leads to a rather more complex evolution, illustrated 
in Fig. [2TJ The variation on the ascending red-giant branch is very similar to the IMq 
evolution illustrated in Fig. 1181 With the helium ignition, the star moves through a 
somewhat convoluted path, shown dotted in the figure, towards larger AIIi, before the 
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Figure 21. 'Hertzsprung-Russell diagram', in terms of the asymptotic dipolar g-mode spacing 
AIIi (cf. Eq. I6.4|l and luminosity, for the 2.5 Mq evolution track illustrated in Fig. [20l The 
dotted part of the curve extends from helium ignition, at a model age of 0.421 Gyr, to the 
establishment of stable helium burning, at 0.439 Gyr. (The end of the track is at 0.565 Gyr.) 



Figure 22. Cyclic buoyancy frequencies N/2-ir in Model 1 (solid). Model 2 (dashed) and 
Model 3 (dot-dashed) in the 2.5 Mq evolution sequence illustrated in Figs 1201 and 1211 

star settles down to quiet helium burning^ During the central helium burning AIIi 
undergoes a further jump to smaller larger values, which is also reflected in the evolution 
track in Fig. [20l 

Again, the variation in AIIi can be understood by considering the evolution of the 
buoyancy frequency, illustrated in Fig. [22l Comparing Models 1 and 2, at roughly the 
same luminosity, the former model is obviously similar to the red-giant IMq model in 
Fig. [191 In Model 2, on the other hand, helium burning causes a convective core where 

f The additional variations along this path appear to be related to modest oscillations in the 
central properties, similar to, but of smaller amplitude, than the oscillations seen in a full-fledged 
helium flash (|SereneUi fc Weiss. 2005 ). 
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the g modes are excluded; since the behaviour of N in the rest of the core is quite similar 
to that of Model 1, the effect is to decrease the integral in Eq. (|6.4|) and hence to increase 
AIIi. This is the dominant cause of the remarkable difference in the period spacing 



betwe en stars on the red-giant branch and in the red clump observed bviBe dding et al 



(|2011l ) and discussed further by Bedding in this volume. The subsequent jump in AIIi 
and L is caused by a jump in the size of the convective core, as illustrated in Fig. 1221 bv 
Model 3 , probably as sociated with changes in the composition profile in the core (e.g., 
[Castellani et al.lll97lh . 



Thi s analysis provides a simple explanation for the behaviour found bv iBedding et al 



(|201lh and furthermore indicates how this may allow the period spacing to be used as a 
sensitive diagnostics of the buoyancy frequency in the stellar core. Further investigations 
will be required to obtain a deeper understanding of the diagnostic power of such data, 
including the ability to probe effects of convective overshoot and rotationally induced 
mixing in the core. 

I am very grateful to Pere Palle and his colleagues at the lAC for the organization of 
an excellent Winter School. I look forward to the next opportunity to participate in such 
a school. J. De Ridder and M.-A. Dupret are thanked for providing several key figures. 
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